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Isotopic effects in the fragmentation of excited target residues following collisions of ^'^C on 
ii2,i24gj^ at incident energies of 300 and 600 MeV per nucleon were studied with the INDRA in 
detector. The measured yield ratios for light particles and fragments with atomic number Z < 5 
obey the exponential law of isotopic scaling. The deduced scaling parameters decrease strongly with 
increasing centrality to values smaller than 50% of those obtained for the peripheral event groups. 
Symmetry term coefficients, deduced from these data within the statistical description of isotopic 
scaling, are near 7 = 25 MeV for peripheral and 7 < 15 MeV for central collisions. 



PACS numbers: 25.70.Mn, 25.70.Pq, 24.10.Pa 

The growing interest in isospin effects in nuclear re- 
actions is motivated by an increasing awareness of the 
importance of the symmetry term in the nuclear equa- 
tion of state, in particular for astrophysical applications. 
Supernova simulations or neutron star models require in- 
puts for the nuclear equation of state at extreme values 
of density and asymmetry Q, |^ . The demonstration 
in the laboratory of the effects of the symmetry term 
at abnormal densities is, therefore, an essential first step 
within a program aiming at gaining such information ex- 
perimentally 

Multifragmentation is generally considered a low- 
density phenomenon, with a high degree of thermaliza- 
tion believed to be reached. Accepting the concept of a 
freeze-out volume and the applicability of grand canon- 
ical logic, the probability of producing a cluster of a 
given atomic number Z and mass A at temperature T 
depends exponentially on the free energy of that clus- 
ter, F{Z,A,T). The cluster free energies depend on the 
strength of the symmetry term Egym = j{A — 2Z)'^/A 
in the liquid-drop energy which, in turn, must depend 
on the extent of expansion of the fragments. This work 



makes use of an observable that isolates the symmetry 
contribution to the cluster free energy to explore the dif- 
ference in this term for fragments produced in periph- 
eral and central collisions. It is found that, while the 
sequential decay strongly degrades the quality of this ob- 
servable, the symmetry energy coefficient does indeed de- 
crease as the collisions producing the fragments become 
more violent. 

In the Copenhagen statistical multifragmentation 
model (SMM), standard coefficients 7 = 23 to 25 MeV are 
used to describe the nascent fragments 0, Qi 13 • 
freeze-out scenario adopted there, normal-density frag- 
ments are considered to be statistically distributed within 
an expanded volume, and the density is only low on av- 
erage. An experimental value for 7 of about standard 
magnitude has recently been obtained within a statisti- 
cal description of isotopic scaling in light-ion (p, d, a) 
induced reactions at relativistic energies of up to 15 GeV 
0. This result, however, may not be representative for 
multi-fragment decays because the data were inclusive 
and the mean multiplicities of intermediate-mass frag- 
ments correspondingly small In the present work. 
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we apply the same method to exclusive data obtained 
with heavier projectiles, ^^C on ^^^Sn and ^^''Sn targets 
at 300 and 600 MeV per nucleon incident energy. Here, 
according to the established systematics [lo| . maximum 
fragment production occurs at central impact parame- 
ters. 

Isotopic scaling, also termed isoscaling, has been shown 
to be a phenomenon common to many different types of 
heavy ion reactions 0jO,E)0|- It is observed by com- 
paring product yields from otherwise identical reactions 
with isotopically different projectiles or targets, and it 
is constituted by an exponential dependence of the mea- 
sured yield ratios R2i{N,Z) on the neutron number N 
and proton number Z of the considered product. The 
scaling expression 

R2i{N,Z) = Y2{N,Z)/Yi{N, Z) ^ C ■ exp{aN + f3Z) 

(1) 

describes rather well the measured ratios over a wide 
range of complex particles and light fragments 0| . 

In the grand-canonical approximation, assuming that 
the temperature T is about the same, the scaling pa- 
rameters a and /3 are proportional to the differences of 
the neutron and proton chemical potentials for the two 
systems, a = A/Xn/T and f3 = A/Xp/T. Of particular in- 
terest is their connection with the symmetry term coeffi- 
cient. It has been obtained from the statistical interpre- 
tation of isoscaling within the SMM and Expanding- 
Emitting-Sourcc Model 1141 and confirmed by an analysis 
of reaction dynamics [l5j . The relation is 

Z^ 

aT = A/ln = /in,2 - Mn.l ~ '^^i^ ~ ^) (2) 

where Zi and Ai are the charges and mass numbers of the 
two systems (the indices 1 and 2 denote the neutron poor 
and neutron rich system, respectively). With the knowl- 
edge of the temperature and the isotopic compositions, 
the coefficient 7 of the symmetry term can be obtained 
from isoscaling. 

The data were obtained with the INDRA multidetec- 
tor in experiments performed at the GSI. Beams of 
^^C with 300 and 600 MeV per nucleon incident energy, 
delivered by the heavy-ion synchrotron SIS, were directed 
onto enriched targets of "^Sn (98.9%) and ^^^Sn (99.9%) 
with areal densities between 1.0 and 1.2 mg/cm^. Light 
charged particles and fragments {Z < 5) were detected 
and isotopically identified with the calibration telescopes 
of rings 10 to 17 of the INDRA detector which cover 
the range of polar angles 45° < 0iab < 176°. These 
telescopes consist of pairs of an 80- /xm Si detector and 
a 2-mm Si(Li) detector which are mounted between the 
ionization chamber and the CsI(Tl) crystal of one of the 
modules of a ring (16i] . Further experimental details may 
be found in and the references given therein. For 
impact-parameter selection, the charged-particle multi- 
plicity Mq measured with the full detector was used, and 
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FIG. 1: Energy spectra of H, He, Li isotopes, and of Be 
and B elements, measured with the calibration telescopes for 
i2q 124 j^gY pgj. nucleon. The dashed hnes indi- 

cate the upper limits set for Z = 1, 2 in the analysis and the 
identification thresholds for Z = 4,5. 



TABLE I: Parameters obtained from fitting the measured 
isotopic yield ratios with the scaling function given in Eq. 







300 MeV 




600 MgV 
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0.0 - 


0.2 


0.28 ± 0.01 


-0.33 ± 0.03 






0.2 - 


0.4 


0.31 ± 0.01 


-0.32 ± 0.01 


0.25 ± 0.02 


-0.28 ± 0.04 


0.4 - 


0.6 


0.36 ± 0.01 


-0.39 ± 0.02 


0.32 ± 0.02 


-0.34 ± 0.03 


0.6 - 


1.0 


0.62 ± 0.01 


-0.68 ± 0.02 


0.52 ± 0.02 


-0.59 ± 0.03 



four bins were chosen for the sorting of the data. For 
600 MeV per nucleon, the two most central bins were 
combined for reasons of counting statistics. 

Kinetic energy spectra of light reaction products with 
Z < 5, integrated over the impact parameter and the 
angular range 6'iab > 45°, are shown in Fig.^ To reduce 
preequilibrium contributions, upper limits of 20 MeV and 
70 MeV were set for hydrogen and helium isotopes, re- 
spectively, which, however, are not crucial. The spec- 
tra of Li, Be, and B fragments were integrated above 
the energy thresholds for isotopic identification which 
amounted to 28, 40, and 52 MeV, respectively. 

The ratios of the fragment yields measured for the two 
reactions and integrated over the chosen intervals of en- 
ergy and angle ( 6'iab > 45°) obey the law of isoscaling. 
This is illustrated in Fig. |21 which shows the scaled iso- 
topic ratios S{N) = i?2i(iV, Z)/exp(/3Z). Their slope 
parameters change considerably with impact parameter, 
extending from a = 0.62 to values as low as a = 0.25 
for the most central event group at 600 MeV per nucleon 
(Table d and Fig.El top). 

Temperature estimates were obtained from the yields 
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FIG. 2: Scaled isotopic ratios S{N) for ^^C + "^'^^*Sn at 
E/A ^ 300 MeV (left panel) and 600 MeV (right panel) for 
the intervals of reduced impact parameter specified in TableQ 
with "central" indicating b/bmax < 0.4 and with offset factors 
of multiples of three. The dashed lines are the results of 
exponential fits according to Eq. Q. Only statistical errors 
are displayed. 



FIG. 3: Isoscaling coefficient a (top), double-isotope temper- 
atures TneLi (middle) and resulting 7app (bottom) for E/A = 
300 MeV (full symbols) and 600 MeV (open symbols), as a 
function of the centrality parameter 1 — 6/6max. The temper- 
atures for the ^^^Sn and ^^*Sn targets are distinguished by 
full and dashed lines, respectively. 



of ^'^He and ^'''Li isotopes, and the deduced ThoLi con- 
tains a correction factor 1.2 for the effects of sequential 
decay 18, 19j. The temperatures are quite similar for the 
two target cases and increase with centrality from about 
6 MeV to 9 MeV (Fig 13 middle). This is consistent with 
the results obtained for ^^^Au fragmentations [3 
and with the established dependence on the system mass 
[ioj. The rise of ThgLI, however, does not compensate for 
the decrease of a, as it did in the case of light-particle 
induced reactions 7], and A/i,i, consequently, decreases 
toward the central collisions. 

The analytical expression for A/in (Eq- contains 
the isotopic compositions of the sources, more precisely 
the difference of the squared Z/A values, A(Z^/A^) = 
(Zi/Ai)^ — (^2/^2)^. For the target spectators, this 
quantity is not expected to deviate significantly from its 
original value , in contrast to mean- field dominated re- 
action systems at intermediate energies fl5. 21]. This was 
confirmed with calculations performed with the Liege- 
cascade-percolation model '22I and with the Relativis- 
tic Mean Field Model (RBUU, Ref. ||23|]) in which iso- 
topic effects of the nuclear mean field are treated ex- 
plicitly. The individual Z/A values change slightly but 
A{Z'^/A'^) remains nearly the same. For central colli- 
sions at 600 MeV per nucleon, e.g., the RBUU calcu- 
lations predict a reduction of A(Z^/A^) by 6% for the 
target-rapidity region after 90 fm/c collision time. 

The suggested corrections are small and can thus be 
temporarily ignored. With the compositions of the orig- 
inal targets, A(Z^/A^) = 0.0367, the expression 7 — 



ar/0.147 is obtained from Eq. [5] It was used to de- 
termine an apparent symmetry-term coefficient 7app, i.e. 
without sequential decay corrections for a, from the data 
shown in Fig. |31 (the mean values were used for T) . The 
results are close to the normal-density coefficient for pe- 
ripheral collisions but drop to lower values at the more 
central impact parameters (Fig.|2| bottom). 

The effects of sequential decay were studied with 
the microcanonical Markov-chain version of the Statis- 
tical Multifragmentation Model (2^]. The target nuclei 
ii2,i24gj^ with excitation energies of 4, 6, and 8 MeV per 
nucleon were chosen as inputs, and the symmetry term 
7 was varied between 4 and 25 MeV. The isoscaling co- 
efficient a was determined from the calculated fragment 
yields before (hot fragments) and after (cold fragments) 
the sequential decay stage of the calculations for which 
standard values for the fragment masses were used. The 
energy balance at freeze-out and during the secondary 
deexcitation was taken into account as described in 3 • 

The hot fragments exhibit the linear relation of a with 
7 as expected (Fig.01 top panel). With 7 = 25 MeV, the 
sequential processes cause a slight broadening of the iso- 
topic distributions and the resulting a is lowered by 10% 
to 20%, similar to what was reported in 0. For smaller 
values of 7, however, larger changes of a are predicted. 
The decay of the wings of the wider distributions of hot 
fragments toward the valley of stability causes the result- 
ing cold distributions to be narrower and the isoscaling 
coefficients to be larger. The overall variation of a with 
7 is weaker, and the decrease of 7 with centrality should, 
thus, be stronger than that of 7app (Fig. 13). For central 
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FIG. 4: Top panel: isoscaling coefficient a for hot (open cir- 
cles) and cold fragments (dots) as a function of the symmetry 
term coefficient 7 as predicted by the Markov-chain calcula- 
tions for "2.124 

Bottom panel: region in the A(ZV^^ )-versus-7 plane consis- 
tent with the measured value a = 0.29 for central collisions 
and with the Markov-chain predictions for cold fragments 
(shaded). The dashed lines indicate the I^{Z'^ jA^) = 0.0367 
of "^'^^''Sn and the RBUU prediction. 



collisions, this implies 7 < 15 MeV as an upper limit 
but much smaller values are more likely to be expected. 
With the present calculations, for excitation energies up 
to 8 MeV per nucleon and the original target composi- 
tions, the measured a < 0.3 is only reproduced with 7 « 

6 MeV (Fig.H top). 

The bottom panel of Fig. 0] shows how the situation 
changes when the variation of the isotopic compositions 
of the two systems is again included as a degree of free- 
dom. The shaded band in the plane of A{Z'^/A^) versus 7 
represents the region consistent with the weighted mean 
a = 0.29 measured for the central bins {b/b„^ax < 0.4) 
at the two energies. It was obtained from the predic- 
tions 0(7) for cold fragments at excitation energies be- 
tween 6 and 8 MeV per nucleon (Fig. 01 top) accord- 
ing to A{Z^/A^) = 0.0367 • 0.29/a(7), i.e. using that 
a is, to first order, proportional to the difference of the 
isotopic compositions, here expressed by A{Z'^ /A^) (see, 
e-g-, If this difference remains close to the cas- 

cade and RBUU predictions the resulting symmetry term 
coefficient for the central reaction channels is very small, 

7 < 10 MeV. To restore the compatibility with a stan- 
dard 7 25 MeV would require considerable isotopic 
asymmetries in the initial reaction phase, much larger 
than what is expected according to the models. 

In conclusion, the observed decrease of the isoscaling 
parameter a with centrality which is not compensated 
by a correspondingly increasing temperature requires a 



decreasing symmetry term coefRcient 7 in a statistical 
description of the fragmentation process. The effect is 
enhanced if sequential fragment decay is taken into ac- 
count. Values less than 15 MeV, as obtained from the 
present analysis, are not necessarily unreasonable in a 
realistic description of the chemical freeze-out state. Be- 
sides the global expansion of the system also the possi- 
bly expanded or deformed structure of the forming frag- 
ments as well as their interaction with other fragments 
and with the surrounding nucleon gas will have to be 
considered. The presented results depend, crucially, on 
the isotopic evolution of the multi-fragmenting system 
as it approaches the chemical freeze-out and, more quan- 
titatively, on the treatment of sequential decay in the 
analysis. These questions deserve further attention. 
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